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Abstract 

Thermal leptogenesis requires the reheating temperature Tr ^ 3 x 10 9 GeV, 
which contradicts a recently obtained constraint on the reheating temperature, 
7r~10 6 GeV, for the gravitino mass of 100 GeV— 10 TeV. This stringent con- 
straint comes from the fact that the hadronic decays of gravitinos destroy very 
efficiently light elements produced by the Big-Bang nucleosynthesis. However, it 
is not applicable if the gravitino is the lightest supersymmetric particle (LSP). We 
show that this solution to the gravitino problem works for the case where the next 
LSP is a scalar charged lepton or a scalar neutrino. We point out that there is an 
upper bound on the gluino ^ 1.8 TeV so that the energy density of 

gravitino does not exceed the observed dark matter density f^DM^ 2 — 0.11. 



1 Introduction 



In the light of experimental data of neutrino oscillations the leptogenesis pQ is the most 
interesting and fruitful mechanism for explaining the baryon-number asymmetry in the 
universe. A detailed analysis on the thermal leptogenesis PJ requires the reheating tem- 
perature Tr ^ 3 x 10 9 GeV, which, however, leads to overproduction of unstable gravitinos 
jH] • Namely, decays of gravitinos produced after inflation destroy the success of Big Bang 
nucleosynthesis (BBN) (HE]. This problem is not solved even if one raises the gravitino 
mass m 3 / 2 up to 30 TeV j^J. Thus, the thermal leptogenesis seems to have a problem with 
the gravity mediation model of supersymmetry (SUSY) breaking. It is, however, pointed 
out in [7] that this gravitino problem may be solved if the gravitino is the lightest SUSY 
particle (LSP). 1 

We show, in this letter, that there is an upper bound on gluino mass, m g i U i no ^6 1.8 TeV, 
for the above solution to work. This is because the heavier gluino produces more abun- 
dantly the gravitino after the inflation and the mass density of the produced gravitino 
LSP, Q 3 / 2 h 2 , exceeds the observed energy density of dark matter, Q^m/i 2 ~ 0.11 
The above bound on the gluino mass will be tested in the next generation accelerator 
experiments such as LHC. 

We show that nonthermal gravitino production by decay of the next LSP (NLSP) 

plays a crucial role of determining the upper bound on gluino mass. We find that a 

consistent NLSP is a scalar charged lepton or a scalar neutrino. Its mass is stringently 

constrained as m 3 / 2 < mNLSP < "Vumo, where the gravitino mass is in the region of 

777,3/2 — 10 — 800 GeV. Other candidates for the NLSP, that is, wino, Higgsino, scalar 

quark (squark) and gluino, are excluded by a new strong constraint from the BBN jlOj. 

except for light gluino of mass m gluino ^ 70 GeV. However, such a light gluino seems to be 

excluded by the present accelerator experiments [TTJ H2] . 

The gauge mediation model provides a solution to the gravitino problem 0. In this letter we consider 
only the gravity mediation model of SUSY breaking. 



2 



2 Cosmological constraints and candidates for NLSP 

As we have mentioned in the introduction, even if the gravitino is the stable LSP, it does 
not totally solve the cosmological gravitino problem in the thermal leptogenesis. Namely, 
we have to constrain its relic abundance consistent with the recent WMAP result jH], 
n 3/2 h 2 < n DM h 2 = 0.11261 

aoi8i- The production rate of the gravitino in the thermal 
bath depends on the reheating temperature T R and the mass-squared ratio ^gi U i no / m 3/2- 
The resultant relic density of the gravitino is calculated as [T31 H3j 
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where a 3 (/i) is a gauge coupling constant of SU(3) C at the scale /i ~ 1 TeV and the gluino 
mass m g i U i no is the one given at the weak scale. 2 The above WMAP constraint on VL^^h 2 < 
Odm^ 2 gives an upper bound on the gluino mass at a given reheating temperature T R and 
a given gravitino mass 1713/2- 

Another important constraint comes from late-time decay of the NLSP into a gravitino 
and its superpartner. The decay width of the NLSP is approximately given by IT1] 



NLSP 



1 m NLSP ^) 
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where m^LSP is the mass of the NLSP and M* = 2.4 x 10 18 GeV the reduced Planck scale. 
In terms of the lifetime, it is written as 
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fi / ^3/2 V /3UUGeV\ 

and then the NLSP can decay during or even after the BBN releasing a large amount 
of energy. Therefore, we have to consider seriously effects of its decay on the BBN to 
confirm the validity of a given model. 

In order not to spoil the success of the BBN, we have to satisfy two constraints on the 
abundance of the NLSP before its decay. One of them comes from the hadronic energy 
release associated with the NLSP decay. Recently, a detailed analysis on the hadronic 



2 Other SUSY particles such as wino and bino contribute to the gravitino production. Thus, the result 
Eq. is regarded as the theoretical lower bound on the gravitino abundance. 



3 



effects has been carried out for a wide range of the NLSP lifetime According to this 
research, even if we take a very conservative bound, the NLSP abundance is constrained 

as 

B h x m NLSP >NLSP < l(T 13 GeV for r NLSP > 10 3 sec . (4) 

Here, is the branching ratio of the NLSP decay into the hadroninc components. Inlsp 
is defined as the yield of the NLSP before its decay, and is given by Ynlsp = 
where n 7 and n NL sp are the densities of the photon and the NLSP, respectively. Further- 
more, if we take the constraint from 6 Li/ 7 Li into account, this constraint becomes much 
stronger as 

^xm NL s P r N LSP<10- {15 ~ 16) GeV for 10 3 sec < r NLSP < 10 8 sec . (5) 

Another constraint comes from the photo-dissociation of light elements by the NLSP 
decay. A detailed investigation on this effect was done in Ref. jH], and the result is 
available in Fig. 2 of Ref. : 

B em x m NLS pF N LSP < KT 12 GeV for r NLSP > 10 6 sec , (6) 

where B em is the branching ratio into electromagnetic components 3 . 

The SUSY standard model (SSM) has various candidates for the NLSP, that is, wino, 
Higgsino, squark, gluino, bino and scalar lepton (slepton). All candidates besides the bino 
and the scalar lepton have dominant hadronic decay modes and hence they are subject 
to the new constraint Eq. (jlj) from the BBN. Our numerical calculation shows that the 



yields of those particles are 

m wino F wino ~ 10- 11 - 2 Gev( T ^^) 2 , (7) 

raffiggsmoYffigssino ~ 10~ 10 ' 7 GeV ( Yo^off^ ) ' ( 8 ) 

m squark F squark ~ 10_11 - 4GeV (T^v) 2 ' ( 9 ) 

m gluino r glumo ~ 10- 12 - 3 Gev(^|^) 2 . (10) 



3 For 10 4 sec^; tnlsp 10 6 sec, this constraint is rather weak; B em x mNLSP^NLSP ~ 10 ^ GeV 0. 
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For our purpose, we have calculated the relic densities of the candidates for the NLSP 
by using micrOMEGAs computer code which includes all possible co-annihilation 
effects. 4 We see that non of them satisfies the cosmological constraint Eq. (j3J except 
for the light gluino of mass ^ 70 GeV. However, as pointed out in the introduction, this 
interesting possibility was already excluded by the present accelerator experiments [TB]. 5 
We should note here that if lifetimes of the NLSP's are shorter than 10 2 sec the 
constraint from the hadronic effects on the BBN becomes weaker ^U] as 

B h x m NL spF NLS p < 4 x l(T 9 GeV for KT 1 sec < t NLS p < 10 2 sec . (11) 

We find that this weaker constraint can be satisfied when the NLSP mass is smaller 
than a few TeV (see Eq. ©-(HH). However, we see from Eq. (JTJ) 6 and © that £lf /2 h 2 
is always larger than Odm'j 2 for tnlsp ~ 10 2 sec and Tr^3 x 10 9 GeV. As a result, no 
NLSP candidate which have dominant hadronic decay modes satisfies all the constraints, 
Eq. JTU), £lf /2 h 2 < tt DM h 2 and T R > 3 x 10 9 GeV, simultaneously. 

On the other hand, the bino and the scalar lepton are still possible candidates for 
the NLSP, since hadronic energy releases from their decays are very small ^3 E| and 
they are not subject to the strong constraint Eq. (0J). However, the bino NLSP is not 
interesting, since its electromagnetic decay violates the constraint Eq. © as pointed out 
in Therefore, we concentrate ourselves to the case of the scalar lepton NLSP in the 
subsequent sections. 

3 Upper bound on the gluino mass with the scalar 
charged lepton NLSP 

There are two candidates for the scalar lepton NLSP. One is the scalar charged lepton and 
the other the scalar neutrino. The scalar neutrino decays into a gravitino and a neutrino. 

4 The above numerical expressions for the NLSP abundance depend on parameters of models and may 
be enhanced by about a factor of 3, which does not, however, affect the following discussion. Here, we 
have assumed that the annihilation processes of the NLSP's do not take place near poles of some particles. 
We have neglected also the nonperturbative QCD effects for the gluino annihilation process. 

5 OPAL and CDF data exclude the existence of (quasi)stable gluino in the mass range 3 GeV 
mgi U i n o 23 GeV and 35 GeV ^ m g i U i no i5 130 GeV 16 . Thus, the mass of the stable gluino is still 
allowed between 23 GeV and 35 GeV. 

6 Hcre, we set m g i u i no — ttinlsp m Eq- for a conservative estimation. 



5 



The effects on BBN from the produced high-energy neutrino will be discussed in the next 
section. 

We consider, in this section, the consequence of the constraints we have discussed in 
the previous section, taking the lightest scalar charged lepton (probably stau) to be the 
NLSP. 7 Constraints for the stau NLSP come from the photo-dissociation of light elements 
by the stau decay, and then we determine the upper bound on the gluino maSS TTigluino 
from the relic gravitino abundance in the following procedure. (In the slepton NLSP case, 
the hadronic contribution dominantly comes from three- and four-body decay channels, 
such as IZG and IqqG ^7j. The branching ratios for these modes are highly suppressed 
as Bh = 10~ 3 -10~ 5 . However, even in this case, if we take the constraint from 6 Li/ 7 Li 
very seriously, the allowed regions that we will present in the following are likely to be 
reduced. In the rest of the paper, we take a conservative point of view assuming this 
constraint to be avoided.) 

Since the relic abundance of the gravitino fl^ 2 h 2 must be smaller than ^dm^- 2 — 0.11, 
we obtain the upper bound on the gluino mass at a given gravitino mass from Eq. ([TJ. 
We show the resultant upper bound on mgi uino in Fig. HJ-(a). We see that this constraint 
does not provide a correct upper bound on m g i U i no , since in addition to the thermal relic 
abundance of the gravitino Q^ 2 h 2 , there is a nonthermal contribution from the late-time 
decay of the stau NLSP, Vtfjfh 2 . 

To determine tt™ 2 T h 2 at a given gravitino mass, we calculate the abundance of the 

stau NLSP, m stau Y st au, before its decay at a given stau mass m stau . Our numerical result 8 

is m stau F s tau ~ 10~ ia3 (m stau /100GeV) 2 for tan/3 = 30, where tan/5 is a ratio of vacuum 

expectation values of the two neutral Higgs bosons (H®, H%), tan/3 = (H®) / (H®). 9 Then, 

7 The possibility that the stau is the NLSP was considered in the framework of gauge-mediation 
model P2|. 

8 Here, we have assumed the gravity-mediated supersymmetry breaking at the grand unified theory 
(GUT) scale ~ 2 x 10 16 GeV. We have set masses squared of the sleptons m 2 lepton + (200GeV) 2 for the 
first and second generations and m 2 lopton for the third generation at the GUT scale, where m s i op ton is a 
universal slepton mass. This boundary condition is favorable to avoid accidental co-annihilations. In this 
sense our estimation on Ystau is conservative. 

9 The above numerical expression for rngtau^stau may depend on parameters of models, although the fol- 
lowing discussion will not change very much. For example, we find m s ta U istau ~ 10 _101 (m st au/100GeV) 2 
for tan/3 = 10. It also may be reduced by a factor of three for the boundary condition which gives a 
larger mass to the right hand scalar lepton than the left one at the GUT scale. 
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we can determine the upper bound on the lifetime of the stau NLSP from the Fig. 2 
of Ref. [5] for a given m atau , which determines the upper bound on the gravitino mass 
^3/2 (see Eq. (JjJJ). By reversing this argument, we can obtain the lower bound on the 
stau mass and hence its abundance, Q stan h 2 , at a given gravitino mass. By converting 
this lower bound on Q stan h 2 to the VL^j^-h 2 by VL^j^-h 2 = (m 3 /2/m stan )Q stan h 2 , we obtain 
the lower bound on the Vfi^j^h 2 a t a given gravitino mass. Finally, we obtain the upper 
bound on the gluino mass for each set of (to3/2,Tr) in order to make the total gravitino 
relic density tt 3 / 2 h 2 = Vtf /2 h 2 + Vt^fh 2 not to exceed the WMAP result, DM /i 2 ~ 0.11. 

The result of the above procedure is given in Fig.[I]-(b), which shows the upper bound 
on the gluino mass for a given gravitino mass. We see that the upper bound reaches 10 
1.3 TeV at m 3 / 2 ~ 200 GeV for the reheating temperature Tr^3 x 10 9 GeV and tan/3 
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Figure 1: The upper bound on the gluino mass at a given gravitino mass. The solid lines 
show the upper bounds on the gluino mass for the reheating temperature Tr = 10 10 GeV, 
3 x 10 9 GeV and 10 9 GeV from the bottom up, respectively. The dashed line denotes the 
lower limit on the stau NLSP mass for a given gravitino mass. We include nonthermal 
relic abundance of the gravitino in the panel (b) (tan/3 = 30). 



Here, we comment on the falling-off behavior of the upper bound in the region of 
m 3/2 ~ 200 GeV. This behavior comes from the fact that the nonthermal production of the 
gravitino from the stau decay becomes dominant. Namely, in the region of 7713/2 ^ 200 GeV 

10 If one adopts the reheating temperature Tr ^ 10 10 GeV for the leptogenesis 0, we find the upper 
bound on the gluino mass to be 600 GeV. 

11 We also find that the upper bound on m g i U i n o reaches 1.1 TeV at m 3 / 2 — 160 GeV for Tr ^ 3 x 10 9 GeV 
and tan/3=10. 
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the gravitino produced in the stau decay dominates over the relic gravitino produced just 
after the inflation. In this sense the mass density of the gravitino is given by the low- 
energy parameters and it is almost independent of the reheating temperature Tr. Thus, 
we find a gravitino DM scenario becomes more attractive in this falling-off region. 

Finally, we should note that the shaded regions in the Fig. ^ are conservative ones, 
thus, they are not always allowed for a given value of (m stau , Mgiuino), since we use the lower 
bound on m stau at a given gravitino mass to estimate Qjw£ T /i 2 . Alternatively, we can obtain 
the upper bound on m g i uino for a given m sta u in the following procedure. As discussed 
above, we can obtain the VL stan h 2 and the upper bound on m 3 / 2 at a given m stau . Then, we 
search the upper bound on m g i uino which satisfies &f/ 2 h 2 + VL^j^h 2 ^ f2 DM /i 2 ~ 0.11 for 
each given m stau within the above gravitino mass bound. We show in Fig. |21 the allowed 
parameter region in the (m s t au , Mgiuino) plane. 
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Figure 2: The upper bound on the gluino mass at a given stau mass (tan/? = 30). The 
solid lines show the upper bounds on the gluino mass for the reheating temperature 
Tr = 10 10 GeV, 3 x 10 9 GeV and 10 9 GeV from the bottom up, respectively. The dashed 
line denotes the lower limit on the gluino mass at a given stau mass, m g i uino = m stau . 



4 Upper bound on the gluino mass with the scalar 
neutrino NLSP 

In this section, we consider the constraints on another candidate for the NLSP, a scalar 
neutrino. A constraint on the scalar neutrino NLSP comes from the destruction of the 
light elements produced at the BBN epoch, which is caused by the high energy neutrino 
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injection from the scalar neutrino decay. The high energy neutrino scatters off the back- 
ground neutrino and produce a lepton-antilepton pair, then it produces many soft photons 
through electromagnetic cascade processes and destructs the light elements. A detailed 
analysis on the effects of the high energy neutrino injection was made in Ref. ^Sj for the 
case of the scalar neutrino LSP. We convert the constraint on the reheating temperature 
in the Fig. 2 in ^H] to the NLSP abundance as 

m NL sp x Fnlsp < l(T 8 GeV for 10 3 sec < r NLS p < 10 9 sec. (12) 

Our numerical result for the abundance of the scalar neutrino is m^Yc ~ 10 _10 ' 8 (mj>/100GeV) 2 , 
where subscript v denotes the scalar neutrino. We see that the scalar neutrino of mass 
^5 3 TeV satisfies the cosmological constraint Eq. (fT2"|) . As we will see in the following 
discussion, this constraint is not significant to determine the upper bound on the gluino 
mass which is given by the requirement for the total gravitino relic density not to exceed 
the WMAP result, tt DM h 2 ~ 0.11. 

To obtain the upper bound on the gluino mass, let us remember that the lower bound 
on the scalar neutrino mass is fixed (by the definition of the NLSP) as m^, > m 3 / 2 at a 
given gravitino mass. 12 Then, it determines the lower bound on the relic abundance of 
the scalar neutrino, Qc,h 2 , at a give gravitino mass. As in the previous section, this lower 
bound on Q„h 2 give rise to lower bound of the Q nonT h 2 at a given gravitino mass. Finally, 
we obtain the upper bound on the gluino mass for each set of (717,3/2, Tr) in order to make 
the total gravitino relic density Q 3 / 2 h 2 = Vtf /2 h 2 + Q^J^h 2 not to exceed the WMAP 
result, fi DM /i 2 ^ 0.11. 

The result of the above procedure is given in Fig. El which shows the upper bound on 

the gluino mass at a given gravitino mass (panel (a)) and at a given scalar neutrino mass 

(panel (b)). We see that the upper bound reaches 1.8 TeV at m 3 / 2 — 450 GeV for the 

reheating temperature Tr ^ 3 x 10 9 GeV. 

12 For the scalar neutrino of mo ^ TO3/2 (the scalar neutrino LSP), the constraint from the direct de- 
tection experiment for the dark matter is very stringent [20J . 



9 




' " 200 400 600 800 200 400 600 800 100012001400 

m 3/2/GeV m v ~/GeV 



Figure 3: The upper bound on the gluino mass at a given gravitino mass (panel (a)) and 
at a given scalar neutrino mass (panel (b)). The solid lines show the upper bounds on 
the gluino mass for the reheating temperature Tr = 10 10 GeV, 3 x 10 9 GeV and 10 9 GeV 
from the bottom up, respectively. The dashed line in the panel (a) denotes the lower limit 
on the scalar neutrino NLSP mass (mo = WI3/2), and the one in the panel (b) denotes the 
lower limit on the gluino mass (m g i uino = mo)- Here, we use the relic abundance of the 
scalar neutrino as m^Y^ ~ 10 _10 ' 8 (m i >/100GeV) 2 . 

5 Conclusions 

A recent detailed analysis jTU] of the hadronic effects of the gravitino decay on the 
BBN leads to a stringent constraint on the reheating temperature of the inflation as 
Tr ^5 10 6 GeV. This constraint contradicts the condition for the thermal leptogenesis, 
that is Tr ^ 3 x 10 9 GeV. A solution to this serious problem is provided p] if the grav- 
itino is stable LSP. However, it depends on nature of the NLSP if this solution works 
or not. We have shown in this letter that the consistent candidate for the NLSP is a 
scalar charged lepton or a scalar neutrino. We have found that there are upper bounds 
on the gluino mass ^ 1.3 TeV and 1.8 TeV, for the former and the latter case, 13 respec- 
tively, so that the density of the gravitino do not exceed the observed dark matter density 
n DM h 2 ~ 0.11 for T fi >3x 10 9 GeV. 

In the present analysis we have used the perturbative calculation in evaluating the 
yield of relic gluino. However, it should be kept in mind that the nonperturbative QCD 

13 We see from Fig. |21 and Fig. that we can not require the universality for three gaugino masses 
at the GUT scale for a large part of the allowed region, since the bino must be heavier than the scalar 
lepton. 
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dynamics may increase the annihilation cross section of gluino which decreases the yield 
of gluino NLSP For instance, if it increases the annihilation cross section by a factor 
of 100, the gluino of mass ^ 300 GeV satisfies the BBN constraint Eq. (4). In this case 
the reheating temperature Tr can be easily taken above Tr ~ 3 x 10 9 GeV. 
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